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It is a hard and important problem to find the criterion of the set of positive-definite matrixes 
which can be written as reduced density operators of a multi-partite quantum state. This problem 
is closely related to the study of many-body quantum entanglement which is one of the focuses 
of current quantum information theory. We give several results on the necessary compatibility 
relations between a set of reduced density matrixes, including: (i) compatibility conditions for 
the one-party reduced density matrixes of any Na X Nb dimensional bi-partite mixed quantum 
state, (ii) compatibility conditions for the one-party and two-party reduced density matrixes of any 
Na x Nb x Nc dimensional tri-partite mixed quantum state, and (iii) compatibility conditions for 

<^^> the one-party reduced matrixes of any Af-partite pure quantum state with the dimension N® M . 

O ! PACS number (s): 03.67.-a, 03.65. Ta 
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I. INTRODUCTION 

Any quantum states (pure or mixed) in finite dimensional Hilbert space can be represented by positive-definite 
matrixes. If the Hilbert space C has a tensor product structure C = C\ ® C2 €5 • • ■ <£> Cm, given a known global density 
matrix P12—M in the space C for the composite system, it is straightforward to calculate the reduced density matrixes 
p a (a = 1, 2, • • • , M) for each subsystem a. However, the reverse of this problem becomes much more involved. Given 
I ' some positive-definite matrixes p a , it is hard to determine whether they can be written as reductions of some global 
iy-) , quantum state. This is what the compatibility problem concerns about. We can list, for instance, the following 
H comparability problems: 

m : 

• Find out the criteria for the one-party density matrixes p a (a = 1, 2, • ■ • , M) so that they can be written as 
reduced density matrixes of some global quantum state p\2---M in the whole Hilbert space C. In this case, we 
typically put some restrictions on the global state P12—M1 for instance, we may require p\i---M to be pure or to 
have a known spectrum (eigenvalues). Otherwise, the trivial assignment P12—M = Pi ® pi ® • ■ • ® PM makes the 
above question meaningless. Pure states actually correspond to a special class of the known-spectrum states 
P12—M with the eigenvalues A{ = 1 and A2 = A3 = • • • = (Aj (i = 1, 2, • • •) are arranged in the decreasing 
order) 

& : 

• Find out the criteria for the two-party (or multi-party) density matrixes p a p (a, (3 = 1, 2, • • • , M) so that they 
can be expressed as reduced density matrix of some global quantum state pn-u- This problem is much more 
involved than the first one, and also much more important for various applications for which we will mention 
several examples. In this case, even without any requirement on the global state pi2---M, in general there is no 
trivial solution to the above problem. 



The investigation of the compatibility conditions between the local and the global quantum states has some impor- 
tant implications: Firstly, this problem is closely related to the study of multi-partite entanglement, which is one of 
the focuses of current quantum information theory [1-6] . The complexity of the compatibility conditions really comes 
from the many-body entanglement inherent in the global quantum state px%--M- If we restrict the global state Px2 - m 
to be separable states, the compatibility condition would be very simple. On the other hand, some general under- 
standing of the compatibility conditions will shed new light on the properties of multi-partite entanglement. Although 
we have achieved remarkable understanding of bi-partite entanglement through consideration of local operations and 
classical communications, we still know little about multi-partite entanglement. Secondly and more importantly, the 
understanding of the compatibility conditions will have significant applications in computational many-body physics 
[8]. In general, the physical interactions are sufficiently local, so the interaction energy typically only depends on some 
reduced local density matrixes. If we find out some general compatibility relations between the reduced and the global 
quantum states, we could design some powerful variational approaches to solve various many-body problems. Finally, 
the study of the compatibility conditions would also help us to get a better understanding of the basic structure of 
quantum mechanics and the associated Hilbert space [7,10,11]. 

What do we know now about the compatibility relations between the local and the global quantum states? There 
are several interesting results concerning the compatibility conditions in some special cases. First, the well-known 
GHJW theorem in quantum information theory [1,9] can be considered as a compatibility condition with the global 
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state restricted to be a pure bi-partite quantum state in the Hilbert space ® , where N denotes the dimension 
of each subsystem. The compatibility criterion in this case is that the two reduced density matrixes should have the 
same spectrum (eigenvalues). There are some generalizations of this result. In particular, the compatibility criteria 
between the one-party reduced density matrixes have been found in some recent works [12,13] if the global quantum 
state is restricted to be a pure state in the Hilbert spaces (C 2 )® M , C 2 ® C 2 ® C 4 , or C 3 ® C 3 ® C 3 . 

In this paper, we provide several new results concerning the compatibility conditions in some more general cases. 
We derive pretty strong necessary conditions for compatibility between the one-party reduced density matrixes for 
the cases that the global quantum state is a bi-partite mixed state in the Hilbert space C Na ® C Nb (the dimensions 
Na and Nb are arbitrary) with a known eigenvalue spectrum { \f B } , or an M-party pure state in the Hilbert space 
(C N )® M . Note that the previously known compatibility conditions for the Hilbert spaces (C 2 )® M and C 2 ® C 2 ® C 4 
(a pure state in C 2 ® C 2 ® C 4 corresponds to a mixed state in C 2 ® C 2 ) can be considered as special cases of the 
results derived here. We also consider the compatibility conditions between the two-party density matrixes for the 
first time, and derive some necessary compatibility conditions when the global state is a general tri-partite mixed 
state in the Hilbert space C Na ® C Nb ® C Nc (the dimensions Na, Nb, and Nc are arbitrary) with a know eigenvalue 
spectrum {Xf BC }. All these results are derived from a unified mathematical method. 

The paper is arranged as follows: in Sec. II, after mention of some mathematical results which are critical for our 
derivation, we derive the necessary compatibility conditions for the reduced density matrices from mixed bi-partite 
quantum states; then in Sec III, we generalize this result and derive the necessary compatibility conditions for the 
one-party and two-party reduced density matrices from general mixed tri-partite quantum states; and finally, in Sec. 
IV, we derive the necessary compatibility conditions for the one-party reduced density matrices from any M-partite 
pure quantum state with the dimension N® M . 



II. COMPATIBILITY RELATIONS OF REDUCED DENSITY MATRIXES FROM MIXED BI-PARTITE 

QUANTUM STATES 

Our methods for derivation of the compatibility relations extensively use two mathematical lemmas. First let us 
summarize these two lemmas: 

In matrix analysis, there is an important theorem which connects the minimization of a matrix with the matrix' 
eigenvalues. That is the content of the following lemma [14] 

Lemma 1: Let A denote an n x n Hermitian matrix, and U denote any n x r matrix with U*U = I r (1 < r ^ n, 
and U* is the adjoint of the matrix U). Then the minimization 

mm(tr(U*AU)) = \\(A) + A|(A) + • • • + Aj(A), (1') 

where the eigenvalues A|(j4), A|(j4), • • • , ^l(A) are arranged in the increasing order. 

In some cases, not all the column vectors of the matrix U are orthogonal to each other. We can get another more 
convenient lemma. 

Lemma 2: Let A denote an n x n Hermitian matrix and U S M„ xs , each column of the matrix U are normalized. 
The columns of U can be divided into two groups, and the columns in the same group are orthogonal each other. 
Suppose the linear dependent number >c = . |it°*w^| is an integer (where uf and are the columns in the first 
and the second group, respectively). Then we have the following conclusion 

min(tr(U*AU)) > £ \}(A) + g \}(A) (1) 

i=l i=l 

where \\ (A), A!(A), • • • , A^(A) arc the eigenvalues of the matrix A and arranged in increasing order. 

Proof: Without loss of generality, we suppose that the vectors u\, tig, • • ■ , uf are the columns in the group I and 
the vectors u\, u\, ■ ■ ■ , u b s _ [ are the columns in the group II . We suppose the eigenvectors of the matrix A are 

Vi , V2 , ■ ■ ■ , v n , corresponding to \\ (A) , A| ( A) , ■ ■ ■ , A^ ( A) , respectively. So all of the vectors in the group I and group 
II can be expanded by the eigenvectors, that is 

ul = auvi + ai 2 v 2 H h ai n v n , 

v-2 = a 2 \v 1 + a 22 v 2 H h a 2n v n , 

uf = ai t ivi + m a v 2 H h ai t „v n . 
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and 



u\ = 0uVi + P12V2 H h /?i„u„, 

"2 = #21^1 + /?22^2 H h /?2nV„, 



= Ps-LlVx + (3 S -U2V2 H V Ps-l, n V n . 

Since the vectors in the same group are normalized and orthogonal to each other, the indexes satisfy I < n and 
s — I < n. Wc can add some vectors into each group to make this group a complete basis of the space. So the 
coefficients must satisfy the requirements Y^j=i \ a ji\ 2 ^ 1 an d Si=i < 1 • Now we can get the formula 



s b 



mm(tr(U*AU)) = min(^^uf Au{) 



i=l j =a 



mm 

u 



EE 1 

j=i i=i 



i=l 



Si=i IAj| 2 )^}(^-) smaller, we must make the coefficients before 



To make the former function X^/=iEi=i \ a vi ~r L^i=i 
the smaller eigenvalues more bigger. Since the constraints Y^j=i l a j'i| 2 ^ 1 an d Si=i l/%| 2 ^ 1j the coefficients before 
the eigenvalues are not more than 2. However, there is another constraint between the coefficients and fy, that is, 

x = X)i j | u f* u 5| 2 • ^ means that the number of the coefficients 2 before the eigenvalues is not more than x. Under 
these constraints, we can get the lower bound of this function, that is, 



mm(tr(U*AU)) > ^Aj(A) + ]T \}(A). 



This is the end of the proof. 

It is need emphasis that s — x is equal to the least number of the linear independent columns when the parameter x 
is a integer. When the number x is not a integer, the similar results can be found easily. We can find at the following 
that the integer situation is enough for our propose. 

The situation that the matrix U can be divided into more than two groups (the columns in the same group are 
orthogonal to each other) is more difficult to deal with. But, for principle, we can get the similar results by carefully 
calculating the linear dependent number x. Though the general case is very complex, We will give a special case in 
the Sec IV.. 

Now we use these lemmas to a multi-partite density matrix to get some relations between the reduced density 
matrix and the global matrix. When the density matrix is a two particle density matrix, and each particle are qubit, 
the original density matrix of the two qubit can be written as 



Pab 



0100,00 ^00,01 ^00, 10 0100,11 

0101,00 aoi.oi aoi.io 0-01,11 

0110,00 dio,oi 0110,10 0110,11 

0111,00 Qin,oi Qin,io 0111,11 . 



We can get the reduced density matrix of pa as 



PA 



0100,00 0100,10 

0110,00 0-10,10 

A a + A 1 



+ 



Ql01,01 0101,11 

an, 01 011,11 



For this case we have the following theorem 

Theorem 1. For the two-qubit density matrix pab, the eigenvalues between pa,Pb an d Pab have the following 
relations 



\\{A) > \\(AB) + XI(AB) 



(2.1) 



Al(B) > \\{AB) + \\{AB) 



(2.2) 
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\\ (A) + x\ (B) > 2X\ (AB) + X\ (AB) + X\ (AB). 



(2.3) 



where A{(A) and X\(B) are the smaller eigenvalue of pA and ps , respectively; X\(AB), X\(AB), X\(AB), X\(AB) 
are the eigenvalues of pab, and they are arranged in increasing order. 

Proof: With the lemma 1, we can get the smaller eigenvalue of the local density matrix pA 



X\(A) = mm[tr(U*p A U)] 

= mm[tr(U*(A Q + A 1 )U)} 

= mm[tr(U*A Q U) + tr(£/*Ai£/)] 



mm[tr([ul 



«oo,oo ioo,io 




Ml 


110,00 aioao 




u 2 _ 



min[ir([ u\ u 2 



)+ir(K <} 



100,00 $ aoo.io $ 




" Ml " 


$ $ $ $ 







110,00 $ «10.10 $ 






. $ $ $ $. 




.0 



^01,01 101,11 




Ml 


111, 01 111,11 




U2 _ 



+tr([ u{ m* ] 



S ioi,oi $ ioi.il 
$ $ $ $ 
$ in,oi $ an. li 





■o - 




Ml 









. 12 . 



where $ means that arbitrary number will make the equality hold, u* means the conjugate of m^. We notice that the 
position of the elements aijj m (i, j, 1,171 — 0, 1) are the same as the position they are in matrix pab- So we choose the 
proper numbers to make the middle matrix is just equal to the density matrix pab- Then 



X[(A) = min [tr[ u\ u* 2 0]p A B 

Ui 



Ml 






min [tr 



M^ u* 2 

M* M* 



PAB 



+ tr([0 m* u* }pab 



Ml 

Ml 

M 2 

M 2 





Ml 



M2 



> min [tr 



u ll u 12 u 13 M 14 

* * * * 

"21 u 22 u 23 u 24 



PAB 



Un U 2 1 

U12 M22 

U13 M23 

L M14 M 2 4 



= X\(AB) + Xl(AB), 



where vectors {«n, M12, M13, M14} and {m 2 i, M2 2 , M23, M24} are orthogonal.. 

With the same reason, we can obtain X\(B) > X\(AB) + X\(AB). It is more important that we can use this method 
to calculate the relations between the smaller eigenvalues of matrix pa and ps, that is 

X\(A) + Xl(B) 



= mmitrWpAUi)] + mm[tr(U; p B U 2 )} 

Ui U2 



= min [tr 



M ll u 12 

u* n 



'12 



PAB 



Mn 

Mn 

M12 

M12 



+ min [tr 

u 2 



u* 21 u* 22 
u* 21 u* 22 



PAB 



m 2 i 

M22 

M 2 i 

m 2 2 



Since the variables in U\ are independent on the variables in U2,we can combine these two matrixes into one matrix 
U, that is 



X\(A) + X\(B) = mm[tr(U* PAB U)} 
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where 



U 



un U21 

Mil U 2 2 

U12 U21 

U12 u 2 2 



(I) 



Now using the lemma 2, we need to calculate the maximal linear dependent number h between unitary Ui and U2. 
Here the linear dependent number k is a constant 1. On the other hand, we can find that this matrix at least has 
three linear independent columns. Then wc get the relation 

\\(A) + X\(B) > 2X\(AB) + Aj(AB) + Aj(AB). 

QED 

These conditions can be viewed as the necessary conditions for the problem whether the single-qubit reduced density 
matrices are compatible with a two qubit density matrix with the eigenvalues {Aj(AB), A|(AB), Xl(AB), X\(AB)}. 
These conditions have already received by Bravyi [13] in another way. Unfortunately, these conditions are not sufficient, 
and the sufficient conditions need another condition 

X\(A) - X\(B) < min{X\{AB) - X\ (AB), X\(AB) - X\{AB)}. 

This condition can not find by our method easily. Now we turn to consider the general case of the two particle 
situation. 

Suppose the dimensions of particle A and particle B are L and N, respectively. Let {Aj(A), X\(A), • • • , A^(A)}, 
{X\(B), Aj(-B), • • • , A] V (B)} and {X\(AB), X\{AB), ■ ■ ■ , X\ N (AB)} be the eigenvalues of the density matrix pa,Pb 
andpAB 1 respectively, and they are arranged in increasing order. Before giving the following theorem, we need 
define majorization relation between two vectors. Let x = {x\, x\, • • • , x\} and y — {y\ , y\, • • • , y^} are n-dimensional 
vectors and the elements are arranged in increasing order. Then we call the vector x is majorized by vector y [15], 
denoted by y >- x, if for each k (k = 1, 2, • • • , n) the following inequality hold 



E4>E 



T 
Vl 



and the equality hold when k = n. The majorization relation have already been extensively used in quantum 
information [16]. If we define the vector X A = {X\{A), X\{A), X\{A)}, X B = {X\(B), A|(B), X^iB)}, 
Kb = {Ef=i A}(AB), Ej=jv+i ^]{AB), ■ ■ ■ , J2-= { l-i )N+ i ^B) } and X B AB = {^ =1 X}(AB), E )t L+1 \](AB), 
• ' " ) J2f=(N-i)L+i ^}(AB)}- Using these definition, we can get the theorem for the general bi-partite case as the 
following 

Theorem 2. As the note before, we get the following relations between the eigenvalues of pa,Pb and pab, 

Xi B >- X A ; (3.1) 

\ab >■ Ab; (3.2) 

k I kN+lL-kl kl 

^Aj(^) + ^A](B)> Aj(AB) + ^A](^),fc = l,2,-..,L-l,/ = l,2,..-,AT-l. (3.3) 

i— 1 j — 1 i—1 j — 1 

Before giving the proof, We need some discussions about these conditions. The majorization relations (3.1), (3.2) 
between the one-party reduced density matrix and the bi-partite density matrix are not just hold for this special 
case. We can see in the following, the majorization is a universal relations between the eigenvalues of reduced density 
matrices and the multi-partite density matrix. This can be viewed as one of the reasons why the majorization relations 
play an important role in the quantum information. The relations (3.3) tell us that some equalities in the former 
two relations can not be hold at the same time except for some special situations. This fact can be viewed as the 
correlation between the different reduced density matrices. When L = N = 2, this theorem is reduced to the theorem 
1. 
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Proof Using the similar method as the qubit case, we can get 



i=l 



mm[tr{U*p A U)] 



= mm[tr(U*A U) + tr(U* A X U) + ■■■+ tr(U*A N U)} 
= min[tr(U* PabU\)\ 

where the unitary matrix U G Mlxu, U*U = I G and Ai is equal to (0| pab |0) . Using the same method which 
used to construct the matrix (J), we can get a unitary matrix U\ G M^M-KkN, we divide this matrix into k blocks, 
and each block is a LN x N matrix which has the following form 



V>N+l,i 



u (L-l)N+l,i 



MJV+1,1 



UN+1,: 



u {L-l)N+l,i 



U( L -l)N+l,i J LN 



(II) 



rxN 



where only the elements (pN + q, q)(p = 0, 1, • • • , L — 1; q = 1, 2, ■ • • , N) in this block are nonzero and i denotes the ith 
block. The first columns in the different blocks are orthogonal and normalized. So all of the columns in the matrix 
U\ are orthogonal and normalized, that is 



kN 



£a,to4)>j>1(ab). 



i=l 



As the same reason, 



J2>J(B)=mm{tr(U*p B U)} 

= mm[tr(U*B U) + tr{U*B x U) + ■■■ + tr(U*B L U)} 



= min[tr(J7 2 * PabU 2 )] 

where the unitary matrix U 2 G Mlnxli- We divide the matrix U 2 into I blocks, and each block is a LN x L matrix 
which has the following form 



G 



Uli 
UNi 



Uli 
UNi 



(III) 



Uli 
U2i 



UNi 



LNxL 



where only the elements ((p— 1)N + q,p)(p = 1, ■ • • , L; q = 1, 2, • • • , N) in the ith block are nonzero. The first columns 
in the different blocks are also orthogonal and normalized. So we can get 



IL 



i=i i=i 

We also need to find the relations between the eigenvalues of pA and ps 



k 



i 



2 Aj (A) + 2 A] (B) 



3=1 



minitriU^pAUi)] + min[tr(C/ 2 * p B U 2 )}. 

Ui U2 



As the lemma 2, we need to calculate the linear dependent number between matrix U\ and C/ 2 . Here it is easy to get 
the linear dependent number x , where k — Y^\=i Sj=i \ {^i l^j' 1 )! ( 1S * ne *th c °l umn m the unitary matrix U\ 
and Vj is the jth column in the unitary matrix U2) is equal to kl. So we get the result 



kN+lL-kl 



kl 



^Al(A)+^A](i?)> ]T Al(AB)+^A](AB). 



i=i 



QED 

We can also find that there are at least kN + IL — kl linear independent columns in matrix U\ and J7 2 . For 
convenience, we let fc > I, we find in the following that at most kl columns in unitary U\ are linear dependent on the 
columns in {/ 2 . At first, wc point out that in each block of unitary U\ at most k columns are linear dependent on the 
columns in unitary L^- Suppose the jth column of the ith block in the unitary U\ is linear dependent on the columns 
in unitary U2, then only the first columns of the block in unitary U2 contribute to the first N elements of the unitary 
U\. So we can get 

ajiVi + a ]2 V 2 H h ctjiVi = Wj. 

where the vector Vi is the first column of the ith block in the unitary U 2 and the vector Wj is a LN vector 
{0, • • • , 0, mj, 0, • • • , 0} T (uij is the jth element in this vector and is equal to the first nonzero element of the jth 
column of ith block and T means transpose). Since Vi(i = 1, 2, • • • , I) are orthogonal each other for different i, and 
they span a Z-dimensional space. Then at most I orthogonal Wi can be get from Vi(i = 1, 2, • ■ • , I) by liner combina- 
tion. That is, at most I columns in the ith block of unitary U\ are linear dependent on the columns in the unitary 
U2- Since there are k blocks, then there at most kl columns in the unitary U\ are linear dependent on the columns in 
the unitary [/ 2 . We will show that this situation can be reached by letting the element W(j_i)jv+i,i in the ith block of 
the unitary U\ be equal to 1 and the other elements zeros. So we can also get the conditions 



kN+lL-kl 



kl 



£\T(B) + 5>}(i?)> Yl Aj(AB) + ^A](AB). 
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These conditions also can be viewed as the necessary conditions to the problem whether the one-party reduced 
density matrices are compatible with the multi-partite density matrix. The method we used to find these conditions 
is very simple, and it is dependent on neither the number of the particles nor the dimensions of the particles. During 
the calculation, the most important thing is to get the maximal linear dependent number x between the different 
unitary matrices. 



III. COMPATIBILITY RELATIONS BETWEEN ONE-PARTY AND TWO-PARTY DENSITY 
MATRIXES FROM TRI- PARTITE MIXED QUANTUM STATES 

There is a density matrix Pabc, where the particle A, B and C arc in L— dimension, M— dimension and 
iV-dimension Hilbert space, respectively. Let {\\(AB), X\{AB), • • • , X\ M (AB)}, {\\(BC), X\(BC), • • • , X\ IN (BC)}, 
{X\(B), Xl(B), • • • , X\ { (B)} and {X\(ABC), X\{ABC), • • • , X\ MN (ABC)} be the eigenvalues of the density matrix 
PABi Pbc i Pb and pabc, respectively, and they are arranged in increasing order. For convenience, We define the fol- 
lowing vectors X AB = {X\(AB), Xl(AB), • ■ • , X{ M (AB)}, X BC = {\\(BC), Xl(BC), • • -,X^ MN (BC)}, X B = {X\(B), 
Xl(B), • • • , Xl(B)}, X ABC = {X\(ABC), Xl(ABC), • • • , X\ MN (ABC)} and X B AB = {Ef=i Aj(AB), £?f i+1 ^B), 

""' Yh=L(M-\)+1 A(AB)}, A| c = {J2iLl A( BC )i Yh=N+1^1( BC )' ""' X^(M-l)jV+l \( BC )}> ^ABC = 

{Ef=i \}(ABC), E-^ +1 \}{ABC), • • • , E™M-i)iv + i \]{ABC) }, Afg c = {^f =1 *){ABC), E f= L+ i *){ABC), 
• ' ' > Y^j=(MN-i)L+i Aj(ABC) }. For this situation, we can get the following theorem. 

Theorem 3. Using the notes defined before, we can get the relations between the eigenvalues of psc , Pab, Pb 
and pabc as 

\ A A lc >■ ^ab (4.1) 

*abc >■ Abc (4.2) 

A| c >- Xb (4.3.1) 
X B AB y X B (4.3.2) 

fiL+r [tN+s fj,LN+Nr+Ls-rs ^LN+rs 

Al(AB) + J2 X }( BC ) > X l( ABC ) + tf(ABC),p<M,0<r<L,0<s<N. (4.4) 

i=l j=l k=l 1=1 

The relation (4.1), (4.2), (4.3) make sure the universality of the majorization relations between the eigenvalues of 
the reduced density matrix and the multi-partite density matrix. The conditions (4.1), (4.2) and (4.3) can obtain as 
the same as in the bi-partite case easily. The relations (4.4) are new relations, and they include more informations 
than the relations (3.3). We can see in the following that we can get some further results of this relations. We only 
need proof the conditions of (4.4). 

Proof: Using the method used in the bi-partite case 

R=fj,L+r S=fj,N+s 

^2 X K AB ) + J2 X}( B C)=mm[tr(U^pABcU 1 )]+mm[tr(U;pABcU 2 )]. 
i=i =i Ul U2 

where U\ is a LMN x RN matrix, every N columns can be viewed as a block. So matrix U\ is divided into R blocks. 
The iih block has the form 
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Uu 

Uu 

Uu 

UN+l,i 

UN+l,i 

UN+\,i 

U(LM-l)N+l,i 

u (LM-l)N+l,i 

U(LM-l)N+l,i 

where only the elements ((k — 1)N + j, j) (k = 1, 2, • • • , LM; j = 1, 2, • • • , N) are nonzero in the ith block, the first 
columns in the different blocks are orthogonal and normalized. The form of this matrix is the same as the matrix U\ 
used in the general bi-partite case. But together with the following unitary matrix J7 2 , we can find that this situation 
is different from the situation in the bi-partite case. In the bi-partite case, there is only one position that both the 
columns in the blocks of (7 2 and the columns in U\ are nonzero. But in this case, it is not. This can substantially 
effect the linear dependent number >t between the unitary matrixes XJ\ and J7 2 . 

The unitary matrix t/ 2 & Mlmnxls, we divide this matrix into S blocks, and each block is a LMN x L matrix 
which has the following form 

uis 

U2,i 
UMN,i 

Uu% 
U2,i 

UMN,i 

Uui 
U2,i 

UMN,i 

where only the elements ((p — 1)N + q,p)(p — 1, • • • , L; q = 1, 2, • • • , MN) in the ith block are nonzero and the first 
columns in the different blocks are also orthogonal and normalized. Now we calculate the maximal linear dependent 
number yc between U\ and J7 2 . We need to point out two facts about these two unitary matrices. 

The first, If R < L and S < N, we can find at most RS columns in the unitary matrix J7 2 are linear dependent on 
the columns in the unitary matrix U\ as the same discussion in the bi-partite case. The second, when R = uL and 
S = fiN, all of the columns in the unitary matrix 11% can be linear dependent on the columns in the unitary matrix U\. 
This can be reached by letting the elements ((j - l)MN + iN + l, 1) in the [(i - l)i + j]th (i = 0, 1, 2, • • • , fj, - 1; j = 
1, 2, • ■ • , L) block of unitary U\ be 1 and the other elements zeroes; at same time letting the elements (k, (i — 1)N + j) 
(k = (i- 1)N + 1, (i- 1)N + 2, ---iN- 1) in the [(i-l)N + j]th (i = 1, 2, • • • , a; j = 1, 2, • • • , N) be nonzero and the 
others zeroes. Using these two facts, when R = uL + r, S = uN + s, there are uLN columns in U2 linear dependent 
on fiLN columns in U\. Since the columns in the same unitary matrix are orthogonal to each other, we need only 
consider the rest columns in the unitary matrix U\ and /7 2 - So at most /iLN + rs columns in the unitary f/ 2 are linear 
dependent on the columns in the unitary matrix U\. QED 

When some eigenvalues of the density matrices are zeros, we can get some stronger relations between the reduced 
density matrix and the global density matrix. 

Theorem 4: Suppose rank(pABc) — LMN — Ls, rank(pBc) = MN — s, rank(pAB) = LM — r and rank{ps) = 
M — t, if r and s satisfy the condition Nr < Ls,then 



(IV) 



(V) 
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r - 1 
*<[—] + ! 



(5) 



is hold, where [x] is the maximum integer which is smaller than x. 

The proof of this theorem is a little technical. It is well known that the rank of a N x N Hermitian matrix is equal 
to N — k (where k is the number of the zero eigenvalues of the matrix). So we need only consider the number of 
the zero eigenvalues of matrix ps- We construct some orthogonal vectors which corresponding to the zero eigenvalues 
from some known matrixes. In the following proof, at first, we point out that the columns in the matrix A (defined 
in the following) are linear combinations of the columns in matrix B (defined in the following) . Then we use this fact 
to construct f 2 -^] + 1 orthogonal vectors which corresponding to zero eigenvalues of matrix ps- 

Proof To proof this theorem, we need to define the matrix A,B and C for density matrix pab,Pbc & n d Pb, 
respectively. The form of the matrix A is the same as the matrix U\ used in the proof of the theorem 3 and made 
the same division. The matrix B is the same as the matrix U2 in the proof of the theorem 3. but we do a different 
division. We viewed each s columns as a block. So the matrix B is divided into L blocks. Each block has the form as a 
LMN x s matrix. In the (i + l)th block, only the elements (iMN+j, k)(j — 1, 2, • • • , MN; k = 1, 2, • • • , s) are nonzero. 
We can note that the nonzero elements are not dependent on the index i, that is, the element (iMN + j, k) = Vjk- 
The form of the (i + l)th block is 



ViMN+1,1 
ViMN+2,1 



ViMN+1.2 

ViMN+2,2 



ViMN+MN,l ViMN+MN,2 



ViMN+l,s 
ViMN+2,s 

ViMN+MN,s 



(VI) 



The columns in the same block are normalized and orthogonal to each other. 

The matrix C can be divided into many blocks, each block is a LMN x LN matrix, further more, we can divide 
each block into L sub-blocks. The (i + l)th sub-block has the following form 



WiMNA 



WiMN+N,l 



WiMNA 



WiMN+NA 



WjMN,l 



WiMN+{M-\)NA WiMN+NA 
WiMN +(M — 1)N A 

WiM N +(M — 1)N A 



(VII) 



LMNxN 



In this sub-block, only the elements (iMN + kN + (k = 0, 1, ■ • • , M — 1; j = 0, 1, • • • , N) are nonzero and they 
are not dependent on j. For each block, the value of the nonzero element are not dependent on the index i. The first 
columns of the different blocks are orthogonal and normalized. Using the lemma 1, we can get the relations of the 
eigenvalues 



E 



X](AB) = mmtr[A*p ABC A], 

A 



E 

k=i 



\\{BC) =mmtr{B*p A BcB], 

B 
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t 

E 



min tr[C* pabcC]. 



Since X^fe=i ^l(^C) = J2k=i ^I(^^C) = 0, there is a unitary transformation between the Ls minimal eigenvectors 
and the columns of the matrix B. Because of Yl^—i ^}{AB) — and the density matrix pabc has only the Ls zero 
eigenvalues, we can see that the columns in matrix A must be a linear combination of the Ls minimal eigenvectors. 
So all of the columns in the matrix A must be a linear combinations of the columns of the matrix B. For convenience, 
we divide each block of the matrix A into L sub-blocks, each sub-block is a MN x N matrix, the (i + l)th sub-block 
has the following form 



UiMN,l 



UiMN+N,l 



UiMN,l 



UiMN+NS 



UiMN.l 



u iMN+{M-l)N,l UiMN+N,l 
u iMN+(M-l)NS 



HMN+(M-1)N,1 



(VIII) 



MNxN 



Now we can find that each columns in the (i + l)th sub-block of the matrix A must be the linear combination of 
the columns of the (i + l)th block of the matrix B. And we take out the first columns of all of the sub-block and 
to find how many columns are linear independent. As the following proof, there are at least + 1 columns are 



linear independent. At first, we denote the column of the (i 
all these vectors in a matrix form as 



l)th sub-block in the (j + l)th block as Vy. We write 



V 01 
V 11 



V 02 

v 12 



(L-2),l 



(L-l),l V(L-1),2 



V 0r 

v lr 



v, 



(L-2),r 



(L-l).r 



(IX) 



Lxr 



The different columns in this matrix are orthogonal and normalized. At first, we consider the columns. There are at 
least one linear independent vector in each column, that is, every elements in the same columns are equal to the same 
vector multiply by a scalar. So we can rewrite the matrix (IX) as the following matrix 



Q!i,iVi 0:2,1^2 

&l,L-lV\ 



a r ^V r 



(X) 



where a^.j (i — 1, 2, ■ • • r, j — 1, 2, • • • , L) are complex number. Now we use the orthogonal conditions that the different 
columns are orthogonal each other. That is, 

(on \ aj )(Vi \Vj) =0i,j = l,2,--.,r 

where the vector |<x,) is ai,2, ••• , c*i,l) T - From this equations we know that if the vectors Vi and Vj are linear 

dependent, then the vectors |aj) and \otj) must be orthogonal. Since the vector |<x,) is in L Dimension space, there are 
at most L Vi are linear dependent to the same vector v p and the vectors \ctj) are orthogonal each other. So there are 
at least + 1 linear independent v p . Using the Schmidt method, we can get l 2 ^] + 1 orthogonal and normalized 

vectors Wi(i = 1, 2, • • • , L). All this vectors can be expressed as the linear combination of the columns of the matrix 
B. 

If we let the first column of the (i + l)th sub-block of the matrix C be equal to one of the vectors Wi. Then the 
columns of the (i + l)th sub-block also can be expressed as the linear combinations of the columns of the (i + l)th 
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block of the matrix B. Since the nonzero elements of the matrix B and C are not dependent on the index of i, the 
columns in the block which the sub-block belongs to can be expressed as the linear combination of the columns of the 
matrix B. So the matrix ps must have at least [^j^-] + 1 zero eigenvalues. 
QED 

From the proof of this theorem, we can find that the most important fact is that the columns in the matrix A 
must be the linear combinations of the columns of the matrix B. Since the symmetry of the particles, exchange the 
role of the density matrix of pab and pbc, the former theorem is hold too. This theorem is the character of the 
multi-partite density matrix, the bi-partite situation has no theorem similar as this theorem. This can be viewed as 
a new correlation between the particles. This theorem is the further results of the relations (4.4). The parameters r 
and s must satisfy the constrain of the inequality (4. 4). 

IV. COMPATIBILITY RELATIONS BETWEEN ONE-PARTY DENSITY MATRIXES FROM N-PARTITE 

PURE QUANTUM STATES 

When using our method to the multi-partite density matrix, we can get much more complicated relations between 
the reduced density matrices and the multi-partite density matrix. The majorization relations between them are hold 
and can be gotten easily. But it is very difficult to find the relations which are similar as the relations (3.3) and (4.4), 
since it is not easy to find the maximal linear dependent number. There is a special case where the TV-partite density 
matrix is a pure state and all of the particles are in a M dimension Hilbert space, we can get some simple formulas, 
that is, 

Theorem 6: For a JV-partite pure state, if every particle are in the M-dimensional Hilbert space, the eigenvalues 
of the one-party reduced density matrices satisfy the following relations 

N M—l p p 

E E A !(i)+E A J( fc )^E A !W' p = i,2,-",M-i,fc^/ = i,2,.--,jv, (6) 

j—lj^k,l i—1 i—1 i—1 

where A} (j), Aj (j), ■ ■ ■ , A] w _ 1 (j) are the eigenvalues of the partial density matrix pj and they are arranged in increasing 
order. 

If N = 2, we can get the necessary and sufficient conditions for the single particle partial density matrices compatible 
with the bi-partite pure state. When M = 2, it can also give the necessary and sufficient compatibility conditions 
between the set of one-party reduced density matrices and the TV-partite density matrix [12]. 

Proof We consider the (N — l)-partite density matrix pi23-N-i- From this density matrix, we can get the one-party 
reduced density matrices pi(i = 1, 2, • • • , N — 1). Now we use the method before to get 

N-2M-1 p N-2 

E E A J0') +E A 1(^ - !) = E ™*MU;p 13 ... N - 1 U i )] + min[tr(C^_ 1 pi2... JV -iC^_i)] 

j — 1 i—1 i—1 i—1 

where the unitary matrix Ui (i = 1, 2, • • • , N — 2) can be divided into M—l blocks, each block is a M" _1 x M N ~ 2 
matrix. The unitary matrix Un-i can be divided into p blocks, each block is also a M N ^ x M N ~ 2 matrix. The 
form of the block (the positions of the nonzero elements in the block) is independent of the block in the same unitary 
matrix Ui. The firsts column in different blocks of the same unitary matrix arc orthogonal and normalized. Further 
more, we can divide the blocks into some sub-blocks. For the block in the unitary matrix Ui, we can divide it into 
M 1 ' 1 sub-blocks, the form of the jth sub- block of the fcth block in Ui is the following 
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l (j-l)M N - i +l,k 



l (j-l)M N ~ i +l,k 



l (j-l)M N - i +M N - i - 1 + l,k 



Uj M N-i_ M N-i-l + lk 



u (j-l)M N - i +M N - i ~ 1 + l,k 



Uj M N-i_ M N-i-l + ltk 



u (j-l)M N - i +l,k 



U(j-l)M N - i +M N - i ~ 1 +l,k 



Uj M N-i_ M N-i-i + i jk 



(XI) 



M N - 1 xM N - i ~ 1 



where only the elements ((j - l)M N ~ i + (I - l)M N ~ i - 1 + k, k) (I = 1, 2, • • • , M; k = 1, 2, • • • , M) are nonzero. The 
nonzero elements of the first columns in the different sub-blocks are equal to each other for the same index I. We can 
find there are some self-similar property on the position of the nonzero elements in unitary matrices Ui and f/j+i- In 
this special case, our main task is to find the number of the least linear independent columns in the unitary matrices 
Ui(i = 2, 3, • • • , N — 1). Using the same discussion in the bi-party case, we find that there are at most (M — l)M N ~ 2 
columns in the unitary matrices Ui(i = 2, 3, • • • , N — 1) are linear dependent on the columns in the unitary U\, and 
at most pM N ~ 2 columns in matrix Un-i are linear dependent on the columns in the matrix U\. This can be reached 
by letting the elements ((k — 1)M N ~ 2 + j, = 1,2, ■ • ■ , M N ~ 2 ) be equal to 1 in the fcth block of the unitary matrix 
U\ . We use the same method to consider the columns remaining in each unitary until the remaining columns only in 
the unitary Un-i- Then we sum up all of the number of the linear independent columns in each unitary matrix, we 
get 

(M - l)M N - 2 + (M - l)M N - 3 + h (M — l)M + p = M N - X -M+p. (7) 

Since the ./V-partite state is a pure state, the nonzero eigenvalues of the density matrix is the same as the nonzero 
eigenvalues of the density matrix p\2—N-i as the Schmidt theorem. So the density matrix P12—N-1 has at most 
M nonzero eigenvalues and the first M^ 1 — M eigenvalues are zeroes. The relations between the eigenvalues of 
the one-party density matrices pi(i — 1, 2, ■ • • , N — 1) and the global density matrix pi2---N-i become the relations 
between the eigenvalues of the one- party reduced density matrices pi(i — 1,2, ■ ■ ■ , N — 1) and pn- So the relations 
between the eigenvalues of the one party density matrix pi(i = 1, 2, • • • , N) are 

N-2M-1 p p 

j—l i—l i—l i—l 

This is the end of the proof of the theorem 6. QED 

From this proof, we can find that the nontrivial relations between the one-party reduced density matrices of the 
multi-partite pure state in this situation must include at least (M — 1)(N — 2) + 1 eigenvalues (the same eigenvalues 
are calculated repeatedly) at the left of the inequality. Further more, using our method we can find almost all of the 
linear relations between the eigenvalues of the one-party reduced density matrixes, only need carefully consider the 
number of the orthogonal vectors needed to express all of the columns. We conjecture that the necessary and sufficient 
conditions of the compatibility problem between the one-party density matrices and a pure state can be expressed 
by the linear relations between these eigenvalues of the local one-party density matrix, that is, the necessary and 
sufficient conditions can form a polytope in the eigenvalues space. If this conjecture is true, we can find the necessary 
and sufficient conditions for the compatibility problem by our method. The necessary and sufficient conditions of the 
compatibility between the single qutrit density matrices and a pure state in C 3 ® C 3 ® C 3 is 

\\(A) + \l(A) < \\(B) + \\{B) + \\{C) + \\{C) 
\\(A) + \\{A) < X\(B) + \\{B) + \\{C) + Xl(C) 
\l(A) + \l(A) < X\(B) + XI(B) + XI(C) + xl(C) 
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X[(A) + 2A^(A) < X[(B) + 2X l 2 {B) + X[(C) + 2Xl{C) (8) 
2X\(A) + Xl(A) < X\(B) + 2Xl(B) + 2X\{C) + X\{C) 
2X\ {A) + A3 (A) < X\ (B) + 2X\ (B) + 2A^ (C) + X\ (C) 
2X\{A) + X\{A) < 2X\(B) + Xl(B) + X\(C) + 2X\{C) 

and the conditions permutation A, B and C. These conditions which are found by Higuchi [12] can be found in our 
method more conveniently. But it is very difficult to prove that this conditions are sufficient. This necessary and 
sufficient conditions are linear and this support our conjecture. When the number of the particle increasing, the 
simplexes of this polytope increase rapidly, the proof used by Higuchi in the three qutrit case is not convenient. We 
must need another method to proof the convex property of this set. 



V. SUMMARY 



In this paper we use a theorem of the analysis matrix to give a simple method to find the relations between the 
reduced density matrix and the multi-partite density matrix. We find hat the majorization relations are the universal 
relations between the eigenvalues of reduced density matrices and the multi-partite density matrix. We also give some 
relations of the eigenvalues between the different reduced density matrixes. All of the relations received in this paper 
can be viewed as the necessary conditions of the problem whether the reduced density matrix is compatible with a 
multi-partite density matrix. What is the necessary and sufficient conditions for the compatibility problem between a 
arbitrary set of density matrices and a multi-partite is far from completely solved, even the special problem whether 
a set of one-party reduced density matrixes is compatible with a pure multi-partite state is very difficult. But the 
method used in this paper give us a possible way to solve these problems, especially for the compatibility problem of 
the pure state. In this paper we only analyzed the number of the linear independent vectors in the N-partite case, if 
we analyze more carefully such as the coefficients of the linear combination, we may find more relations. This is one 
of our further work. On the other hand, if we add some symmetry such as the translation invariant on the particles, 
we can get more constrains on the eigenvalues of the density matrices, then we may find the necessary and sufficient 
conditions of the compatibility problem. This symmetry is very useful in condensed matter and the quantum phase 
transition. This problem will also be investigate in the future. 
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